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In this paper, we study the thermal entanglement in a two-qubit Heisenberg XYZ system with different 
Dzyaloshinskii-Moriya (DM) couplings. We show that different DM coupling parameters have different in- 
fluences on the entanglement and the critical temperature. In addition, we find that when Ji (i-component 
spin coupling interaction) is the largest spin coupling coefficient, Di (i-component DM interaction) is the most 
efficient DM control parameter, which can be obtained by adjusting the direction of DM interaction. 
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I. INTRODUCTION 



Entanglement has been studied intensely in recent years due 
to its fascinating nonclassical feature and potential applica- 
tions in quantum information processing QJJ] . As a simple sys- 
tem, Heisenberg model is an ideal candidate for the genera- 
tion and the manipulation of entangled states. Many physical 
systems, such as nuclear spins JjL quantum dots yfl, super- 
conductor |0] and optical lattices J5[], have been simulated by 
this model, and the Heisenberg interaction alone can be used 

. Recently, the 
XY model I 



for quantum computation by suitable codin 
Heisenberg models, including Ising model _ 
XXX model U, XXZ model d and XYZ model HI E 
have been intensively studied. Shan et al. investigated the ef- 
fects of DM interaction, impurity and exchange couplings on 
entanglement in XY spin chain it 1 311 - Ay diner et al. studied the 
thermal entanglement of a two-qutrit Ising system with DM 
interaction Il4tl . they found that the control of entanglement 
can be optimized by utilizing competing effects of the mag- 
netic field and the DM interaction. Wang et al. investigated 
the effects of the DM interaction and intrinsic decoherence on 
entanglement teleportation in the two-qubit XXX Heisenberg 
model J3. 

In the above models, the influences of the z-component 
DM interaction (arising from spin-orbit coupling) and the 
external magnetic field on the entanglement have been dis- 
cussed, but the DM coupling interactions along other direc- 
tions have never been taken into account. Quite recently, 
we discussed the influences of x-component DM interaction 
on entanglement in Heisenberg XXZ model 

d and XYZ 

model 1 17]. To research further the differences between DM 
coupling interactions along different directions, in this pa- 
per, we generalize the special Heisenberg models to the gen- 
eralized Heisenberg XYZ models with different DM inter- 
actions, and then analyze the different influences of D x (x- 
component DM control parameter), D y (y-component DM 
control parameter) and D z (z-component DM control param- 
eter) on the entanglement and the critical temperature. We 
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find that Di is the most efficient DM control parameter when 
Ji is the largest spin coupling coefficient. Thus, according 
to the relation among Ji(i = x,y,z), we can know which 
is the most efficient DM control parameter. In order to pro- 
vide a detailed analytical and numerical analysis, here we take 
concurrence as a measure of entanglement 118!]. The con- 
currence C ranges from to 1, C = and C = 1 indi- 
cate the vanishing entanglement and the maximal entangle- 
ment respectively. For a mixed state p, the concurrence of 
the state is C(p) = max{2A ma x — St=i A,-, 0}, where A^s 
are the positive square roots of the eigenvalues of the matrix 
R = p(a y (g)o- y )p*(a v 0cr y ), and the asterisk denotes the 
complex conjugate. 

This paper is organized as follows. In Sec. II, we intro- 
duce the Heisenberg XYZ models with different DM inter- 
action parameters, and give the analytical expressions of the 
concurrences. In Sec. Ill, we analyze the different influences 
of different DM control parameters (D x , D y and D z ) on the 
entanglement and the critical temperature. Finally, in Sec. IV 
a discussion concludes the paper. 



II. THE HEISENBERG XYZ MODELS WITH DIFFERENT 
DM INTERACTION PARAMETERS 



A. Heisenberg XYZ model with D x 

The Hamiltonian H for a two-qubit anisotropic Heisenberg 
XYZ chain with DM interaction parameter D x is 

H = J x afa 2 ° + J v a\a 2 J + J z a z (T z + D x (a\o- z .-a z a%), (1) 

where Jj(t = x, y, z) are the real coupling coefficients, D x is 
the x-component DM coupling parameter, and a l (i = x,y, z) 
are the Pauli matrices. The coupling constants Ji > corre- 
sponds to the antiferromagnetic case, and Ji < corresponds 
to the ferromagnetic case. This model can be reduced to some 
special Heisenberg models by changing Jj. Parameters Ji and 
D x are dimensionless. 

In the standard basis {|00), |01), 1 10) , |11)}, the Hamilto- 



2 



nian ( 1 ) can be rewritten as 



H = 



( Jz H-^X J X Jy 

~J% Jx ~h Jy ^Dx 
Jx Jy Jz ~^J^x 



V Jx Jy 



-iD, 



Jz 



By calculating, we can obtain the eigenstates of H: 

|*i> = ^(|01> + |10», 



(3a) 



As.4 = -^e" 



-Jx±Jy=fJz 



(6b) 



where Z = 2e~^ cosh(^A) + 2e* cosh(p). Thus, the 
^ concurrence of this system can be expressed as jflal : 

f max{|Ai - A 3 | - A 2 - A 4 ,0}, if J y >J z , m 
S max{|Ai - A 4 | - A 2 - A 3 ,0}, if J y < J z . U> 



which is consistent with the results in Ref. [ 161 for J r = J, 



'v 



|vI/ 2 ) = — (|00) + |11)), 



(3b) 



1 



|* 3 ) = (sin 6»i 1 00) -i cos 6>i 1 01) +« cos (9i 1 10) -sin 6»i 1 1 1) ) , 

(3c) 

|*4) = -^(sin0 2 |OO)-Hcos0 2 |Ol)-icos0 2 |lO)-sin0 2 |ll)), 
v2 

(3d) 

(4a) 
(4b) 



with corresponding eigenvalues: 

E\,2 — Jx ± Jy T Jz, 



E3A = - Jx ± W, 



where 0\ 2 



arctan( 



2D, 



■), and w 



\/4:Dx + (J y + Jz) 2 - The system state at thermal equi- 

exp( t7 H t ) 

librium (thermal state) is p(T) = — , where 

Z = Tr [cxp( j^gT )] * s tne P art iti° n function of the system, 
H is the system Hamiltonian, T is the temperature and 
Kb is the Boltzmann costant which we take equal to 1 for 
simplicity. Thus, in the above standard basis, we can get the 
following analytical expression of the density matrix p(T): 



B. Heisenberg XYZ model with D y 

Here we consider the case of the two-qubit anisotropic 
Heisenberg XYZ chain with y-component DM parameter D y . 
The Hamiltonian is 

H' = J x a¥a%+J y a\a%+J z alaZ+D y (ataZ-<Tf<TZ), (8) 

where D y is the y-component DM coupling parameter, which 
is also dimensionless. 

In the standard basis {|00>, |01), 1 10) , |11)}, the Hamilto- 
nian (8) can be rewritten as 



H' = 



Dqi Jr J II 



I Jz 

D y ~Jz Jx ~l~ Jy 

JJy Jx Jy Jz 

V Jx Jy &v J^ 11 Jz 



-D„ 



(9) 



,y J^y 

Similarly, by calculating, we can obtain the eigenstates of if': 



1*1) 



1 
71 



(|01) + |10)), 



|* 2 ) = -^(|oo)-|n», 



(10a) 



(10b) 



where 



mi 2 



P(T) 



I m 1 q q* m 2 
q* m n 2 q 
q n 2 ni q* 
m 2 q* q m\ 



— (e T ±e T sin 0i ± e T sin 2 J 



1 / -El . 2 n 1 2 n \ 

H,2 = 7T^\ e T ±e t cos 0\ ± e t cos 2 ), 



|*3) = ^(sin0i|OO)-cos^i|Ol)+cos0i|lO)+sin^i|ll)), 
(5) (10c) 

= ^(sin0 2 |OO)-cos(/)2|Ol)+cos02|lO)+sin ( /)2|ll)), 
v2 

(lOd) 
(11a) 
(lib) 



with corresponding eigenvalues: 

^12 = Jy ^ Jx ~f J zi 



E' 3 A = -J y ± W', 



where 6\ 2 = arctanf T 2 "P" — 7), and u>' = 
4Z?2 + ( J x + J z ) 2 . In the above standard basis, the 
density matrix p'(T) has the following form: 



(7= — (e T sin 0i cos 0i — e r sin 2 cos 02 J. 
*2Z 

After straightforward calculations, the positive square roots of 
the eigenvalues of the matrix R = p(a v (g) a y )p*(cr y a y ) 
can be expressed as: 



1 i±. 
Ai, 2 = — e t 



(6a) 



p'(r) 



mi -g' g' m 2 

-g' rai n 2 -g' 

q' n' 2 n'i q' 

m! 2 —q' q' m! Y 



(12) 



3 



where 

/ 1 i , -Si- -Sk • 2 i 



--Si • 2 i -\ 
e sin 02 



1 , _£i 



-^i 2 
- ' e T ± e T cos 



6x ± e ^ cos 2 </>2 ) , 



— (e r sm 0i cos 0i 



Then the positive square roots or me eigenvalues c 
R' = p'(a y (jy)p'*(<jy a v ) can be obtained 



e T sin 02 cos 02). 
roots of the eigenvalues of the matrix 

II O _7/\ 1 1 



1 -J„±JxTJ» 

A ll2 = t 



A 3,4 = 



1 Jy±" 



(13a) 
(13b 



where Z' = 2e-r L cosh(iy^) + 2e^ cosh(^). Thus, the 
concurrence of this system can be expressed as: 

f ma^Ai - A 3 | - A' 2 - X' 4 ,0}, if J x > Jz, nA 
- \ maxllAa-A^I-Ai-A^O}, i/ J x < J,. 114 




FIG. 1 : (Color online) (a) The concurrence is plotted versus D x (blue 
solid line) and D y (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D x — 3 (blue solid 
line) and D y = 3 (red dashed line). Here the coupling constants 
J x = 0.2, J y = -1 and J z = -0.5. 




C. Heisenberg XYZ model with D z 

The Hamiltonian H" of a two-qubit anisotropic Heisenberg 
XYZ chain with z-component DM parameter D z is 

H" = J x olo x 2 + J v a\al + J z a{a z 2 + £> z (<r*crf - crfcr*), 

(15) 

where D z is the z-component DM coupling parameter, which 
is also dimensionless. 

Using the similar process, we can get the eigenstates of H": 



l< 2 ) 



1*3,4) 



1 
1 

V2 



(|00) ± |11)), 
(|01>±x|10», 



with corresponding eigenvalues: 

El 2 — Jz i J x =P t/y, 



where % 



££ 4 = -J z ±w", 

J;c+Jy-2lD z 



(16a) 



(16b) 



(17a) 



(17b) 



and w h 



\/4D| + (J x + Jj,) 2 - Similarly, we can get the analyti- 
cal expressions of p(T)" and R" , but we do not list them 
because of the complexity. After straightforward calcu- 
lations, the positive square roots of the eigenvalues of 
R" = p"{a y a y )p"*(a y a y ) can be expressed as: 



1 Jk±s 

e T 



(18a) 



FIG. 2: (Color online) (a) The concurrence is plotted versus D x (blue 
solid line) and D y (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D x = 3 (blue solid 
line) and D y = 3 (red dashed line). Here the coupling constants 
Jx = -1, J H = 0.2 and J z = -0.5. 



1 

3,4 



(18b) 



with Z" = 2eT^ cosh(^y^) + 2eT cos h(%). Thus, the 
concurrence of this system can be written as: 



f max{|A'/ - A 3 '| - A 2 ' - A£, 0}, i/ J. 

— S r 1 a // \ // 1 \ // \ // /at - p t 



> J, 



max{|Ai' - X'l\ - A 2 ' - A3, 0}, */ J x < J y 



(19) 



which is also consistent with the results in Ref. Il6ll when 

J X Jy. 

From Eqs. (7), (14) and (19), one can see that when J x = 
J y , there is the same entanglement for D x = D y ; when J y = 

and when 
= D z . So 



J z , there is the same entanglement for D y = D z 
Jx = Jz, there is the same entanglement for D x 



when J x = J v 



J z , there is also the same entanglement for 



the same values of DM interaction parameters {D x , D y and 
D z ) 



III. THE COMPARISON BETWEEN THE DIFFERENT DM 
CONTROL PARAMETERS IN HEISENBERG XYZ MODEL 

A. The comparison between D x and D y 

In Heisenberg XYZ model, we has analyzed all kinds of 
spin coupling coefficients satisfying J. x > J y . For simplicity, 
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FIG. 3: (Color online) (a) The concurrence is plotted versus D y (blue 
solid line) and D z (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D y = 3 (blue solid 
line) and D, — 3 (red dashed line). Here the coupling constants 
J x = -0.5, J v = 1 and J z = 0.2. 



(a) (b) 
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FIG. 4: (Color online) (a) The concurrence is plotted versus D y (blue 
solid line) and D z (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D y — 3 (blue solid 
line) and D, — 3 (red dashed line). Here the coupling constants 
J x = -0.5, J y = 0.2 and J z = 1. 



here we choose one set of spin coupling coefficients satisfying 
Jx > Jy, and plot Fig. 1 to demonstrate the properties of dif- 
ferent DM parameters. In Fig. 1(a), we find the entanglement 
increases with the increase of DM coupling parameter. Fur- 
thermore, the critical value of D x is smaller than D y , and D x 
has more entanglement for D x = D y . In Fig. 1(b), we find 
that increasing temperature will decrease the entanglement, 
and D x has a higher critical temperature than the same D y , so 
that the entanglement can exist at higher temperatures for D x . 

Similarly, we has analyzed various spin coupling coeffi- 
cients satisfying J x < J y . For simplicity, we choose one set 
of coupling coefficients satisfying J x < J y , and plot Fig. 2 
to demonstrate the properties of different DM parameters. In 
Fig. 2(a), it is shown that increasing the DM coupling pa- 
rameter can enhance the entanglement. Besides, for a certain 
temperature, the critical value of D y is smaller than D x , and 
D y has more entanglement for D x = D y . In Fig. 2(b), it is 
easy to see that D y has a higher critical temperature than the 
same D x , so that the entanglement can exist at higher temper- 
atures for D y . 

Thus, the x-component parameter D x has a smaller critical 
value, higher critical temperature and more entanglement than 
the same y-component parameter D y for J x > J y , and D y has 
a smaller critical value, higher critical temperature and more 
entanglement than the same D x for J x < J y . 

B. The comparison between D y and D z 

Similarly, for J y > J z case, Fig. 3 is plotted to show 
the properties of different DM parameters in Heisenberg XYZ 
model. In Fig. 3(a), we can see that the y-component parame- 
ter D y has a smaller critical value, and more entanglement for 
D y = D z . In Fig. 3(b), we can see that D y has a higher crit- 
ical temperature than the same D z , so that the entanglement 
can exist at higher temperatures for D y . 

Contrarily, for J y < J z case, the concurrence versus differ- 
ent parameters is shown in Fig. 4. In Fig. 4(a), it is easy to 
see that the z-component parameter D z has a smaller critical 
value, and more entanglement for D y = D z . In Fig. 4(b), it 
is easy to see that D z has a higher critical temperature than 
the same D y , so that the entanglement can also exist at higher 



temperatures for D z . 

Thus, the y-component parameter D y has a smaller critical 
value, higher critical temperature and more entanglement than 
the same z-component parameter D z for J y > J z , and D z has 
a smaller critical value, higher critical temperature and more 
entanglement than the same D y for J y < J z . 

C. The comparison between D x and D z 

Here, for J x > J z case, we plot Fig. 5 to illustrate the prop- 
erties of different DM parameters in Heisenberg XYZ model. 
In Fig. 5(a), D x has a smaller critical value and more entan- 
glement for D x = D z . In Fig. 5(b), D x has a higher critical 
temperature than the same D z , so that the entanglement can 
exist at higher temperatures for D x . 

For J x < J z case, we plot Fig. 6 to illustrate the properties 
of different DM parameters. In Fig. 6(a), for a certain temper- 
ature, D z has a smaller critical value and more entanglement 
for D x = D z . In Fig. 6(b), D z has a higher critical tempera- 
ture than the same D x , so that the entanglement can also exist 
at higher temperatures for D z . 

Thus, the x-component parameter D x has a smaller critical 
value, higher critical temperature and more entanglement than 
the same z-component parameter D z for J x > J z , and D z has 
a smaller critical value, higher critical temperature and more 
entanglement than the same D x for J x < J z . 

The above results indicate that when Ji is the largest spin 
coupling coefficient, the i-component DM control parameter 
Di has the smallest critical value, the highest critical temper- 
ature and the most entanglement. So according to the relation 
among spin coupling coefficients, we can know which is the 
most efficient DM control parameter. 

IV. DISCUSSION 

We have investigated the thermal entanglement in a two- 
qubit Heisenberg XYZ system with different Dzyaloshinskii- 
Moriya (DM) couplings. We find that different DM inter- 
action parameters (D x , D y and D z ) have different influ- 
ences on the entanglement and the critical temperature. When 
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Ji(i = x, y, z) is the largest spin coupling coefficient, the i- 
component DM interaction Di(i = x, y, z) has the smallest 
critical value, the highest critical temperature and the most 
entanglement. In addition, when J x = J y = J z , there is the 
same entanglement for the same values of DM interaction pa- 
rameters (D x , D y and D z ). Thus, according to the relation 
among spin coupling coefficients (J x , J y and J z ), the most 
efficient DM control parameter can be obtained by adjusting 
the direction of DM interaction. 



FIG. 5: (Color online) (a) The concurrence is plotted versus D x (blue 
solid line) and D z (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D x = 3 (blue solid 
line) and D, = 3 (red dashed line). Here the coupling constants 



Jx 



-0.2, J v = 0.3 and J z 




FIG. 6: (Color online) (a) The concurrence is plotted versus D x (blue 
solid line) and D z (red dashed line) for T = 3. (b) The concurrence 
is plotted as a function of the temperature T for D x = 3 (blue solid 
line) and D z = 3 (red dashed line). Here the coupling constants 



Jx 



-l,J y =0.3 and J z = -0.2. 



Acknowledgments 



This work is supported by National Natural Science Foun- 
dation of China (NSFC) under Grant Nos: 60678022 and 
10704001, the Specialized Research Fund for the Doctoral 
Program of Higher Education under Grant No. 20060357008, 
Anhui Provincial Natural Science Foundation under Grant 
No. 070412060, and the Key Program of the Education De- 
partment of Anhui Province under Grant No. KJ2008A28ZC. 



[1] Phys. World 11, 33 (1998), special issue on quantum informa- 
tion; C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 2881 
(1992); C. H. Bennett etal, ibid. 70, 1895 (1993); A. K. Ekert, 
ibid. 67, 661 (1991); M. Murao et al, Phys. Rev. A 59, 156 
(1999). 

[2] B. E. Kane, Nature (London) 393, 133 (1998). 
[3] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998); G. 
Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 59, 2070 

(1999) ; B. Trauzettel etal., Nature Phys. 3 192 (2007). 

[4] T. Senthil et al, Phys. Rev. B 60 4245 (1999); M. Nishiyama et 

al, Phys. Rev. Lett. 98 047002 (2007). 
[5] A. S0rensen and K. M0lmer, Phys. Rev. Lett. 83, 2274 (1999). 
[6] D. A. Lidar, D. Bacon, and K. B. Whaley, Phys. Rev. Lett. 82, 

4556 (1999); D. P. DiVincenzo etal, Nature (London) 408, 339 

(2000) ; L. F. Santos, Phys. Rev. A 67, 062306 (2003). 

[7] D. Gunlycke et al, Phys. Rev. A 64, 042302 (2001); Z. Yang et 

al, Phys. Rev. Lett. 100, 067203 (2008). 
[8] G. L. Kamta and Anthony F. Starace, Phys. Rev. Lett. 88, 

107901 (2002); X. G. Wang, Phys. Rev. A 66, 034302 (2002); 

Y. Sun et al, Phys. Rev. A 68, 044301 (2003); Z. C. Kao et al, 

Phys. Rev. A 72, 062302 (2005); S. L. Zhu, Phys. Rev. Lett. 96, 

077206 (2006). 



[9] M. Asoudeh and V. Karimipour, Phys. Rev. A 71, 022308 
(2005); G. F. Zhang, Phys. Rev. A 75, 034304 (2007). 

[10] G. F. Zhang and S. S. Li, Phys. Rev. A 72, 034302 (2005); Y. 
Zhou et al, Phys. Rev. A 75, 062304 (2007); M. Kargarian et 
al, Phys. Rev. A 77, 032346 (2008). 

[11] D. C. Li and Z. L. Cao, Optics Communications, in press; L. 
Zhou etal, Phys. Rev. A 68, 024 301 (2003); G. H. Yang et al, 
e-print arXiv:quant-ph/0602051 

[12] F. Kheirandish et al, Phys. Rev. A 77, 042309 (2008); Z. N. 
Gurkan and O. K. Pashaev, e-print arXiv:quant-ph/0705.0679 
and arXiv:quant-ph/0804.0710. 

[13] Chuan-Iia Shan etal, Chin. Phys. Lett. 25, 817 (2008). 

[14] C. Akyiiz, E. Aydiner, and O. E. Miistecaplioglu, Optics Com- 
munications 281, 5271 (2008). 

[15] Liang Qiu, An Min Wang, and Xiao Qiang Su, Physica A 387, 
6686 (2008). 

[16] D. C. Li, X. P. Wang, and Z. L. Cao, J. Phys. Condens. Matter 

20 325229 (2008). 
[17] D. C. Li and Z. L. Cao, Eur. Phys. J. D 50, 207 (2008). 
[18] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997); 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 



